Non-classical rotational inertia in the supersolid state 
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An abrupt drop in the moment of inertia recently found 
in solid helium 4 is explained in terms of dynamics of zero- 
point vacancies (ZPV). Mechanical decoupling of ZPV from 
the motion of the container due to Bose statistics is developed 
to a macroscopic phenomenon by repulsive interaction. It 
gives a negative answer to the question whether BEC is a 
necessary condition for non-classical rotational inertia in a 
bulk three-dimensional system. 

PACS numbers: 67.80.-s, 67.40.-w, 05.30.-d 



Solid helium 4 has been termed a quantum crystal. 
On introduction of defects, it displays several anomalies 
not observed in classical crystals. Vacancies can move 
about unusually rapidly. Hence, it is natural to expect 
that the quantum-mechanical fluctuation delocalizes va- 
cancies at a sufficiently low temperature, which is called 
'zero-point vacancies' (ZPV) [1]. In such a 'supersolid 
state', non-classical rotational properties are expected, 
which was argued with emphasis on the connectivity of 
the wave function [2] . A number of experiments searching 
for anomalies in thermodynamical or mechanical proper- 
ties have been performed, but ended with null results [3]. 
Recently, an abrupt drop in the moment of inertia was 
found in the torsional oscillation measurements on solid 
helium 4 confined in a porous media [4] and on a bulk 
solid helium 4 [5] . This discovery raises the question why 
other measurements until now have found null results, 
and leads us to reconsider the definition of superfluidity 
and that of solids. 

A natural way to discuss superfluidity in a confined sys- 
tem is to focus on its rotational properties. Anomalous 
behaviors of the system are normally attributed to the 
Bose condensate; in other words, off-diagonal long-range 
order (ODLRO). The question whether the Bose conden- 
sate is a necessary condition for superfluidity has been 
discussed in terms of academic interest. For solid helium 
4, however, this question is not an academic but a practi- 
cal one. ZPV is still a hypothetical object, but is highly 
probable if solid helium 4 shows superfluidity [3] . A fun- 
damental feature of crystals lies in their periodicity in 
density; that is, diagonal long-range order (DLRO). One 
has to face a serious question whether crystals remain 
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stable while showing superfluidity that violates their pe- 
riodicity. 

In this paper, we assume the existence of ZPV obey- 
ing Bose statistics, but do not assume its Bose condensate 
[6]. Since vacancies are defects on the crystal lattice, it is 
unlikely that crystals with a macroscopic number of de- 
fects remain stable [7]. Rather, in light of the definition 
of solids, a more natural idea is that the number of ZPV 
is large, but not macroscopically large in a solid helium 
4. Hence, a direction to be explored is the possibility 
that under the influence of Bose statistics a microscopic 
number of ZPV's create the non-classical rotational iner- 
tia without ODLRO. For the dynamical properties such 
as superfluidity, the repulsive interaction generally en- 
hances singular properties due to Bose statistics. When 
the container slowly rotates at low temperature, some mi- 
croscopic regions emerge in solids, in which ZPV and sur- 
rounding atoms decouple from the motion of container, 
thus remaining at rest. The repulsive ZPV's are likely to 
spread uniformly in coordinate space. This feature makes 
ZPV outside of these regions behave similarly with ZPV 
at rest. Otherwise, the density of ZPV would become lo- 
cally high, thus raising the interaction energy. The me- 
chanical behavior of microscopic fractions of solids due 
to Bose statistics (decoupling) is developed to a macro- 
scopic phenomenon by the repulsive interaction. In this 
paper, we stress that the non-classical rotational inertia 
(NCRI) is possible even in the case of no condensate by 
such a mechanism, and propose its formalism [8]. 

Zero-point vacancies. Vacancies are usually regarded 
as localized objects which occasionally move from one 
position to another in crystals. It takes an energy ei for 
an atom to diffuse to the surface. Further, the lattice 
distortion due to the vacancy raises the energy by 62. In 
quantum crystals, however, the delocalization of the va- 
cancy due to the quantum-mechanical fluctuation lowers 
the energy by A. Hence, it takes eq = ei + £2 ^ A to 
make a vacancy in the crystal. 

In quantum crystals, there is a case in which A slightly 
exceeds ei-|-e2, thus leading to eo < 0. The possible states 
of such a vacancy are classified by the quasi momentum p. 
With the tight-binding model applied to ZPV, its energy 
spectrum has a form as e(p) = eo -|-p^/2to*, where m* 
is an effective mass of the vacancy involving surrounding 
atoms in its motion. 

When one vacancy approaches another in crystals, a 
lattice-distortion pattern by one vacancy is normally not 
consistent with that by another vacancy. (An atom can 
not follow two displacement patterns that contradict each 
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other.) This situation results in an effective repulsive 
interaction U between the two vacancies. 

Let us consider a crystal with ZPV a.s H = Hq + Hz 
where Hq represents a hamiltonian of a crystal with no 
vacancies, Hz a hamiltonian of ZPV 
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and $ denotes an annihilation operator of ZPV as a 
spinless boson [9] . The motion of vacancies is expressed 
by Green's function as G{u},p) = 1/{lo — e{p)). Multiple 
scattering between vacancies is expressed by t-matrix as 
t = {7/(1 + UG). In the tight-binding picture of ZPV, 
the band width of ZPV has an order of A, thus G ~ 
A~^. As Z7 — > 00, t approaches 1/G ~ A. Hence, an 
energy e of the repulsive ZPV takes a form such as e = 

— |eo| + nA, where n is a number density of ZPV. The 
increase in the number of ZPV stops when e reaches zero, 
that is, n = |eo|/A, following which the number of ZPV 
may change only through their fusion to the surface of 
crystal. Its relaxation time is much longer than the time 
required to establish equilibrium for a given number of 
ZPV. At low temperature, ZPV behaves as a repulsive 
Bose quasiparticle in the equilibrium state. 

For the stability of solids, even in the case of ZPV 
(eo < 0) A only slightly exceeds ei + £2, which would 
otherwise lead to the collapse of a solid. Because of |eo| = 
|(ei + £2) — A| <C A, the number density satisfies n = 
|eo|/A <^ 1. Hence, it is unlikely that the number of 
ZPV is a macroscopic one. 

Moment of inertia. By the torsional oscillator experi- 
ment, one measures the moment of inertia / of a sample 
as a response to the restoring force of the torsion rod. As 
long as the external force is small, / is determined by the 
rotational property of the system without that force. 

Consider a solid helium 4 in a uniform rotation round 
z-axis with angular velocity fl. A hamiltonian in a co- 
ordinate system rotating with the body is H — ft ■ L, 
where L is the total angular momentum. The rotation 
is equivalent to the application of the probe acting on 
the sample. The perturbation —Q L is cast in the form 

— ^^(J7 xr)-p , where flxr = Vd{r) is the drift velocity 
of the sample at point r. With the current density J{r), 
the perturbation has a form as 
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One assumes the spatial homogeneity of the sample dur- 
ing rotation. Hence, J(r) is approximated as 



J(r) = 
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Using Eq,(3) and Vd = (— Oy, fia;, 0) in the right-hand 
side of Eq.(2), and fl = (0, 0, 0) in its left-hand side, one 
obtains the angular momentum L^, hence the moment of 
inertia Iz = Lz/^ as 



4 = X^(0,0) / {x-'+y^)dv. 
Jv 
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(A) The susceptibility of a normal solid satisfies 
X-^(0,0) — x'^(0,0). Hence, replacing x"^(0, 0) by p in 
Eq.(4), one obtains the definition of classical 

The generalized susceptibility is defined as 

Xnu{q, w) = ^^X^{q, t^) + ( V - ) X^iQ, t^)- (5) 



For the later use, we define a term proportional to q^Qi, 
in x^i^ as x^^ so that 



Xu.'^iQ, = Sf.^x'^iq, uj) + q^q^, 
= S,_cuX^ + XiJ.i^{q,uj), 
and rewrite Eq.(4) as 



X^(9,t^) -X'^(9,t^) 



Iz=I^ - lim 
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(x^ + y^)dr. (7) 



(B) At the onset of ZPV, according to iJ = + 
Hz, the susceptibility xiq-^) splits into xa{q,'-^) and 
X^^^ {q,oj) [11]. The normal part gives no contribution 
to NCRI in Eq.(7) due to Xo(O'O) = Xo (O'O). whereas 
the validity of x^-^^'^(0,0) = x^'^^'^(0,0) must be ex- 
amined. For the appearance of NCRI, the balance be- 
tween the longitudinal and the transverse excitation must 
be destroyed on a macroscopic scale {q — > 0). On this 
point, one can trace back to Feynman's physical argu- 
ment on how the Bose-statistical coherence suppresses 
the transverse excitation, and destroys this balance [12]. 
In his argument, not the Bose condensate but Bose statis- 
tics is essential. 

(C) To see the effect of the repulsive interaction, let us 
consider Xi.iv of ZPV as an ideal Bose gas. Within the 
linear response, it is defined as 



Because of divvd{r) = 0, Eq.(2) says that the rotation 
represented by Vd{r) acts as a transverse- vector probe to 
the excitations in solids [10]. We need the linear response 
of the system J(r) to Vd{r). Since a solid normally ro- 
tates like a rigid body, we can use the mass density p in 
J{r) = pVd{r). Microscopically, however, one must be- 
gin with the transverse susceptibility x'^il: of t^i^ sys- 
tem. (By definition, the mass density p is regarded as the 

longitudinal response to a force as p = nm = x^(0,0).) where Ji^{q,T) is a current of ZPV {h = 1) [13] 



ZPV 
A, HI/ 







1 r" 
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and 15 > is a ground sate of Hq + e(p)^p^'p- Its 
Xfiiy{<l,^) is given by 



E 



/(e(p))-/(e(p + g)) ^ 
OJ + e{p) - e{p + q) 



(10) 



where f{e{p)) is the Bose-Einstein distribution. 

(1) If ZPV would form the condensate, f(e(p)) in 
Eq.(lO) is a macroscopic number for p = and nearly 
zero for p 7^ 0. Thus, in the sum over p in the right-hand 
side of Eq.(lO), only two terms corresponding to p = 
and p = —q remain, with a result that 
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where no is the number density of the condensate. Equa- 
tion. (7) with Eq.(ll) shows NCRI. (The relative degree 
{Iz — II)/ 11 of the deviation from the classical moment 
of inertia would be simply proportional to the number 
density of the condensate.) 

(2) When ZPV forms no condensate, however, the sum 
over p in Eq.(lO) is carried out by replacing it with an 
integral, and one notices that dependence disappears 
in the result. Hence, without the repulsive interaction, 
BEC is a necessary condition of NCRI. 

Non-classical rotational inertia. We will formulate a 
mechanism of NCMI under the more realistic condition 
for a solid helium 4. Instead of \S > in Eq.(8), a ground 
state |G > of Ho + Hz must be used as 



< G\TrJ^{x,T)MO,0)\G> 



< S\TrJ„{x,T)X{0,0)exp 
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where Hi{t) represents the repulsive interaction. 
Due to the repulsive interaction U contained in 
exp{—i f Hi{T)dT), the scattering between ZPV's fre- 
quently occurs as illustrated by an upper bubble with 
a dotted line in Fig. 1(a). The current- current response 
tensor J^(a;, r) Jjy(0, 0) is depicted by a lower bubble. 

Since |G > is a ground state under the influence of 
Bose statistics, the perturbation must be developed in 
such a way that as the order of the expansion increases, 
the susceptibility gradually includes a new effect due to 
Bose statistics. An important feature is that ZPV in 
Jfiix, r) J^(0, 0) and ZPV in the upper bubble due to the 
repulsive interaction form a coherent wave function as a 
whole. Hence, we must seriously consider the influence 
of Bose statistics on the graph like Fig. 1(a). When one 
of the two ZPV's in the lower bubble and that in the up- 
per bubble have the same momentum (p = p'), and the 



other ZPV in the two bubbles have the same momentum 
(p + q = p' + q'), a, graph made by interchanging these 
two types of ZPV's must be included in the expansion of 
X^v^- The interchange of ZPV lines with p and p'{= p) 
in Fig. 1(a) yields Fig. 1(b). Further, the interchange of 
ZPV lines with p + q and p + q'{= p + q) yields Fig. 1(c), 
a graph linked by the repulsive interaction, whose con- 
tribution to x^v^ is given by 



f{e{p))^f{e{p + q)) 



■ e(p) - e(p + q) 



(13) 



With decreasing temperature, the coherent wave func- 
tion grows to a large size, and the interchange of ZPV due 
to Bose statistics like Fig.l occurs many times. Hence, 

one can not ignore the higher-order term x^if which 
is significant in the larger coherent wave function. Among 

many ZPV's with various momentums, ZPV's at rest 
with p = play a key role, leading to the following form 



X^r'^"n9,0) = ^/7"FM9r+S 
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where 



Mq) = 



(exp(-/3/z) - l)-i - (exp(/3[e(g) - ^^]) - l)-^ 



<q) 



(15) 



and (U is a chemical potential. Ff){q) is a positive 

monotonously decreasing function of q^ which approaches 
zero as — > 00. As ^ 0, Fp{q) increases at any q. An 
expansion form of Fi){q) around g^ = is given by 



FM = 



/3exp(/3/x) 



(12) 



/J2 ^ {1 + eMm 

- ^exp(/3/i)- 

2 (1 - exp(/3/i))'^ 
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When q^Xjiv I qiiqv has a finite value at g — > in Eq.(7), 
that is, when a power series in U 

00 

x^r(9,o) = ^Et^"^M'?r'' (17) 

diverges as q~'^ at q 0, the moment of inertia shows 
a non-classical behavior. In normal solids (/3/x <C 0), a 
small Fp[ci) guarantees the convergence of X^v^ {,q,^), 
with a result that, 



X^.""(9,0) 



4 l-UF0{q)- 



(18) 



With decreasing temperature (/i 0), however, a grad- 
ual increase of Ff}[q) makes the higher-order term signifi- 
cant in Eq.(17), finally leading to the divergence of power 
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scries in x^l^^ {q,0). The convergence condition is first 
violated at q ■ 



when UFp^O) = 1, that is, 



(19) 



Prom now, we call Tq satisfying Eq.(19) an onset temper- 
ature of NCRI. Equation (19) implies that a solid with 
ZPV shows NCRI without the condensate (/x ^ 0). At 
T = To, substituting Eq.(16) in Eq.(18) with the aid of 
Eq.(19), we find as a leading term 



x^r(9,o) = 



2sinh|/3o/i| 



In view of Eq.(6), Eq.(20) leads to x 



ZPV.L 



X 



ZPV,T 



(20) 

(<?,0) ^ 

(g, 0) at ^ 0. Since the large-scale interchange 
of bosons plays an essential role for deriving Eq.(20), 
this result confirms Fcynman's intuitive argument on the 
influence of Bose statistics on the excitations. Using 
Eq.(20) in Eq.(7) with the aid of Eq.(19), we obtain 



I = p - 

— j-v 



m 



[{x^+y^)dr, (21) 
Jv 



which implies that {Iz — Iz)m fo^' superfluidity by ZPV 
does not allow the one-body interpretation as Ps/ p, but 
reflects the many-body effect due to Bose statistics and 
the repulsive interaction. Since ZPV is not a particle in 
the literal sense, but represents the motion of surround- 
ing atoms, this is a natural result. 

Discussion. Superfluidity is a complex of phenomena, 
and therefore has some different definitions such as, (1) 
persistent current without friction, (2) the Hess-Fairbank 
effect, (3) quantized circulation, (4) almost no friction on 
moving objects in the system below the critical velocity, 
(5) peculiar collective excitations and (6) the Josephson 
effect [14]. (Strictly speaking, the experiment by Kim 
and Chan adds a new definition of superfiuidity to the 
list.) Since these phenomena (1) ~ (6) are observed to 
occur together in most materials exhibiting superfluidity 
(including superconductors), it is difficult to elucidate 
the exact nature of the relationships between these def- 
initions. From the standpoint of this paper, there is no 
reason to expect that all superfluid-like behaviors corre- 
sponding to these definitions should occur in quantum 
crystals. We must therefore classify them into two cate- 
gories. The first is the phenomenon which an interplay of 
Bose statistics and the repulsive force creates without the 
condensate, and the second is that for which the conden- 
sate is necessary. The first and the second corresponds to 
an observable and an unobservable phenomenon in quan- 
tum crystals, respectively. (Prom this viewpoint, we must 
put a new interpretation on some 'null results' in earlier 
measurements of mechanical properties of a solid helium 
4 [15].) We can regard a solid helium 4, together with 
a helium 4 film [16], as appropriate material for examin- 
ing the relationship between the different definitions of 
superfluidity. 
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FIG. 1. The first-order Feynman diagram of the cur- 
rent-current response tensor. The black and the white small 
circle represents a vector and a scalar vertex, respectively. 
The dotted line denotes the repulsive interaction. 
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